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ABSTRACT 

We present the Fundamental Plane (FP) in the z = 0.28 cluster of galaxies RX J0142. 0+2131. There is no 
evidence for a difference in the slope of the FP when compared with the Coma cluster, although the internal 
scatter is larger. On average, stellar populations in RX J0142. 0+2131 have rest-frame V-band mass-to-light 
ratios (M/Ly) 0.29 ± 0.03 dex lower than in Coma. This is significantly lower than expected for a passively- 
evolving cluster formed at Zf = 2. Lenticular galaxies have lower average M/Ly and a distribution of 
M/Ly with larger scatter than ellipticals. Lower mass-to-light ratios are not due to recent star formation: our 
previous spectroscopic observations of RX J0142. 0+2131 E/S0 galaxies showed no evidence for significant 
star-formation within the past ~ 4 Gyr. However, cluster members have enhanced a-element abundance ratios, 
which may act to decrease M / Ly. The increased scatter in the RX J0142. 0+2131 FP reflects a large scatter 
in M/Ly implying that galaxies have undergone bursts of star formation over a range of epochs. The seven 
easternmost cluster galaxies, including the second brightest member, have M / Ly consistent with passive evo- 
lution and Zf = 2. We speculate that RX J0142. 0+2131 is a cluster-cluster merger where the galaxies to the 
east are yet to fall into the main cluster body or have not experienced star formation as a result of the merger. 
Subject headings: galaxies: clusters: individual: RX J0142. 0+2131 - galaxies: evolution - galaxies: stellar 
content 



1. INTRODUCTION 

The Fundamental Plane (FP) of elliptical (E) and lenticular 
(SO) galaxies is a log-linear relation between effective radius, 
surface brightness, and velocity dispersion {e.g. Djorgovski & 
Davis 1987; Dressier et al. 1987; J0rgensen et al. 1996, here- 
after JFK1996). Its ubiquity in nearby clusters of galaxies has 
been noted by many investigators (e.g. Bender et al. 1992; 
JFK1996; Colless et al. 2001). Coefficients in the observed 
FP differ from those predicted by the virial theorem, and this 
'tilt' implies that mass-to-light ratio (M/L) is not a constant 
for E/SOs but depends to some extent on galaxy properties 
(see Cappellari et al. 2006). The FP must confront studies 
showing that the growth of galaxies in clusters involves com- 
plex processes like galaxy mergers, cluster mergers, and AGN 
activity (e.g. Fabian et al. 2000; Tran et al. 2005). 

If E/S0 galaxies form in a short burst and their stellar pop- 
ulations subsequently evolve passively, one would expect a 
gradual decrease in M/L with cosmic epoch. This induces 
an offset in the FP as a function of redshift, with no change in 
slope. However, the tendency of star formation to occur over a 
longer period in lower-mass galaxies (e.g. Juneau et al. 2005; 
Treu et al. 2005), would be manifest as a steepening of the 
empirical FP with redshift. Evidence for changes in the slope 
of the FP are seen in the latest studies of the FP at z > 0.8 
(e.g. Holden et al. 2005; J0rgensen et al. 2006), but this has 
not been observed at 0.2 < z < 0.8 (e.g. van Dokkum & 
Franx 1996; Kelson et al. 2000; Wuyts et al. 2004), most 
likely because these studies do not probe far enough down 
the mass function. The internal scatter of the FP reported at 
0.2 < z < 0.8 is consistent with that of the Coma cluster 
(JFK1996). Since scatter in the FP is a reflection of the varia- 
tion in M/L, this is taken as evidence that cluster E/SOs at a 
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given mass are coeval. 

In this letter we examine the FP in RX J0142.0+2131, a 
cluster of galaxies at z = 0.28 for differences in slope and 
internal scatter with respect to that of Coma. We use Ho = 70 
km s _1 Mpc -1 , O m = 0.3, Qa = 0.7, and convert previous 
work to these values where necessary. The look-back time to 
RX J0142.0+2131 in this cosmology is 3.2 Gyr. 

2. RX J0142.0+2131: PREVIOUS OBSERVATIONS 

We presented high-signal-to-noise GMOS-N spectroscopy 
of 43 spectroscopic targets in the field of RX J0142.0+2131 
in Barr et al. (2005; hereafter BDJBC). Velocity dispersions 
and line indices for 30 cluster members were derived. We 
determined that scaling relations between metal indices and 
velocity dispersion, and the strengths of the 4000A break, are 
inconsistent with a scenario in which stellar populations form 
at Zf > 2 and evolve passively to look like cluster galaxies 
at low redshift (i.e. without additional star formation and/or 
merging). Stellar populations in RX J0142. 0+2131 have a- 
element abundance ratios ([a/Fe]) which are enhanced by, 
on average, 0.14 ± 0.03 over Coma galaxies. Luminosity- 
weighted mean ages for galaxies in RX J0142. 0+2131 are 
similar to those in Coma, i.e. older than would be expected 
at z = 0.28. The cluster velocity dispersion is much larger 
than its X-ray luminosity or richness suggest. Despite this, 
no sign of substructure was found by testing the distribution 
of spectroscopically-confirmed cluster members in RA,dec,z 
space. 

3. DATA REDUCTION AND ANALYSIS 

The data reduction and analysis method for the GMOS- 
N spectroscopy was presented in BDJBC. In this let- 
ter we restrict the surface photometry calculations to 
spectroscopically-confirmed cluster members. We will 
present effective radii, mean surface brightnesses and mor- 
phological classifications for a larger sample at a later date. 

RX J0142.0+2131 was observed with the Hubble Space 
Telescope Advanced Camera for Surveys (ACS) on UT 2003 
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November 01 and 2004 July 03. Two positions were imaged 
in the F775W band for 4420s each. Reductions are performed 
with the PyRAF task MULTIDRIZZLE using the standard 
procedure. ACS surface photometry is a nalysed using th e 
two-dimensional fitting program GALFIT (Peng et al. 2002). 
We fit galaxies' luminosity profiles as as r 1 / 4 profiles. Back- 
ground or companion galaxies may influence the derived pa- 
rameters, so these are fit simultaneously as Sersic profiles 
(Sersi c~968l) . We also classify galaxies by Hubble type us- 
ing a method analogous to that of Smail et al. (1997). Mag- 
nitudes and surface brightnesses are corrected to rest-frame 

V using F775W magnitudes and GMOS colors in a simi- 
lar way as used in BDJBC. A complete description of our 
data, reduction techniques, analysis of potential systematics 
as well as our method of deriving surface-brightness param- 
eters from the GALFIT fitting program will be presented in 
a future paper (Chiboucas et al., in preparation). Twenty- 
eight spectroscopically-confirmed cluster members were cov- 
ered by the two ACS fields. Of these 11 are E, 12 are SO and 
the remaining 5 are spirals, irregulars, mergers, or not classi- 
fiable. 

We estimate the systematic and random errors in the values 
returned by GALFIT using Monte Carlo simulations. Artifi- 
cial galaxies are generated with GALFIT and with the IRAF 
routine ARTDATA/MKOBJECTS. These objects cover the 
full range of parameter space in magnitude, r c , axis ratio,and 
position angle of the RX J0142. 0+2131 sample. The Sersic 
index is varied randomly between 2.5 and 5 and we adjust 
the diskiness/boxiness of each artificial galaxy. Noise is then 
added to the artificial galaxies which are placed randomly in 
the ACS images and recovered as r 1 / 4 profiles. 

The results of the simulations indicate that logr c derived 
for the r 1 / 4 profile is systematically greater than the simulated 
value by 0.06 dex. The quantity, logr e + /31og(/) c , which 
enters the FP and referred to hereafter as the Fundamental 
Plane Parameter (FPP), is known to be much less sensitive 
to the imposition of the r 1 / 4 -law fit (e.g. Lucey 1997). The 
derived FPP is systematically offset from its simulated value 
by —0.02, with a standard deviation of 0.03. The magnitude 
of the deviation is not a function of the FPP. We do not make 
corrections to the surface photometry but quantify the total 
uncertainty on the FPP as 0.05. 

4. THE COMPARISON SAMPLE 

Our comparison sample consists of 116 E/S0 galaxies in 
the Coma cluster; velocity dispersions are from J0rgensen 
(1999). We use Gunn r surface photometry of E/S0 galax- 
ies fit using an r 1 / 4 law from J0rgensen et al. (1995), as well 
as new B and R c photometry obtained with the McDonald 
Observatory 0.8m telescope. We transform fro m Gunn r to 
the rest-frame V using r — R c = 0.354 (J0rg ensen 19941) . 

V = R c + 0.337(5 - R c ) + 0.089 derived for E/S0 galax- 
ies using Bruzual & Chariot (2003) models, a median Coma 
red-sequence color of B — R c = 1.485, and V rcs t — V — 
Qv — 101og(l + z) where the fc-correction, Qy = 0.043, 
is derived from Bruzual & Chariot (2003) models. The sur- 
face brightness parameter in the V^-band FP is log(/) c = 
-0.4« M ) e - 26.43). 

We recalculate the FP for the Coma cluster by minimising 
the absolute residuals in log r c , log tr, log(/) c perpendicular to 
the plane, following JFK1996. Three emission-line galaxies 
are excluded from the fit. We find 

logr c = (1.29 ± 0.07) logcr - (0.83 ± 0.03) log(/) c - 0.28 



with an rms of 0.08 in log r e . 

5. THE FUNDAMENTAL PLANE IN RX 10142.0+2131 

We fit the FP to E/S0 galaxies in RX J0142.0+2131 in the 
same way as for the Coma cluster. Emission-line galaxies 
identified in BDJBC are excluded. The FP is 

logr c = (0.93 ± 0.26) logcr - (0.97 ± 0.11) log(/) c + 1.10 

with an rms of 0.10 in logr c . 

Selection effects can be crucial when comparing two sam- 
ples at varying epochs. As a test we recalculate the FP in 
Coma and RX J0142.0+2131 for the subsamples with My < 
-19.7 and M > 10 10 - 3 M Q . There are no significant differ- 
ences in the coefficients of log a or log(/) c , or rms between 
either the low- or high-redshift sample and its parent sample. 

In order to assess the significance of any difference in slope 
between the Coma FP and that of RX JO 1 42 . 0+2 1 3 1 , we make 
500 random realisations of a high-z cluster. Twenty-three 
galaxies are taken from the Coma cluster and scattered so that 
therms in log r c is equal to that in RX J0142.0+2131. We then 
fit the FP as described above. Values of the coefficients in the 
FP are consistent with Coma 89% of the time, and consistent 
with RX J0142.0+2131 55% of the time. We conclude, there- 
fore, that there is no strong evidence for a difference in the 
slope of the FP for RX J0142.0+2131 when compared with 
Coma. The galaxies from RX J0142. 0+2131 are plotted on 
the Coma cluster FP in Figure 

Even with the imposition of the same mass cutoff at M > 
1O 1O ' 3 M0, the internal scatter of the points from the FP re- 
mains higher for RX J0142.0+2131 than it is for Coma (it 
actually increases to 0.13 in logr e for RX J0142. 0+2131 as 
opposed to 0.08 in Coma). This phenomenon is not seen in 
previous studies of the FP in clusters at similar redshift (e.g. 
Kelson et al. 2000; Fritz et al. 2005). The scatter in the FP 
reflects the variation of Al/L within stellar populations, and a 
high scatter could imply that there are two or more groups of 
galaxies with different mean values of M/L. In Figure^we 
separate galaxies according to E or SO classification; the scat- 
ter for SOs is higher, 0. 12 in log r c compared with 0. 10 for Es. 
We also plot cluster members east of 01h42m06s including 
the second brightest member as filled boxes; open boxes are 
galaxies in the main body of the cluster. The internal scatter 
is 0.05 and 0.12 in log r c for the filled points and open points 
respectively. 

Observables are converted to mass and MjLy via, M = 
5a 2 r c /G, and log MjLy = 2 log a - logr c - log(/) + T, 
where for a in km s , r c in kpc, and (I) e in L Q pc -2 , T = 
—0.73. Figure|2]shows log MjLy vs. logM. 

We determine the slope of the log MjLy vs. logM di- 
agram for Coma and the complete RX J0142. 0+2131 sam- 
ple by minimising absolute residuals perpendicular to the 
slope. The slopes are consistent, though the uncertainty 
for RX J0142.0+2131 is large. We therefore calculate the 
offset in MjLy ratio for RX J0142.0+2131 by preserving 
the Coma slope and calculating the median offset of the 
RX J0142.0+2131 galaxies. This is -0.29 ± 0.03 dex in 
MjLy for the whole sample. Lenticular galaxies have lower 
average MjLy than ellipticals (0.06 dex), and their scatter 
in M/Ly is 0.18 dex as opposed to 0.14 dex. The offset pre- 
dicted by the single-stellar population (SSP) models of Maras- 
ton (2005) with passive evolution and Zf orm = 2 is —0.10, 
shown as the dot-dashed line on Figure [2] If the differences 
in mass-to-light ratio are caused by different luminosity- 
weighted mean ages of stellar populations, then the Maraston 
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FIG. 1 .— (a) The FP for RX J0142.0+2131 and for Coma seen face on. 
Triangles are E/S0 galaxies in Coma, squares Es in RX J0142. 0+2131, and 
circles are SOs in RX J0142. 0+2131. Filled squares or circles are members of 
RX J01 42. 0+2131 east of 01h42m06s. Error bars show the median error for 
each sample. The dashed lines indicate the selection effect of the magnitude 
limits. Dotted lines denote the region not occupied by luminous ellipticals in 
Bender et al. (1992). (b) The FP seen edge-on against the effective radius, (c) 
The FP projected along one of its shortest edges. In (b) and (c) The solid line 
is the best fit to the low-z sample minimising absolute residuals (see text). 



models suggest that galaxies in RX J0142. 0+2131 incurred 
a major star-formation episode at z = 0.85 ±0.1, making 
them only ~ 4 Gyr old. This is highly unlikely given the 
spectroscopic evidence for old stellar populations presented 
in BDJBC. 

6. THE HISTORY OF E/S0 GALAXIES IN RX J0142.0+2131 

The internal scatter of the RX J0142.0+2131 FP is 
larger than that of Coma. This suggests that galaxies in 
RX JO 142. 0+2131 have undergone bursts of star formation 
over a range in epochs. The fact that the SO galaxies have 
lower Mj Ly and larger scatter is consistent with results sug- 
gesting they form stars at lower redshift, or are more recent 
additions to clusters (e.g. Treu et al. 2003; Mcintosh et al. 
2004). 
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FIG. 2. — Mass-to-light ratios inferred for galaxies in RX J0142.0+2131 
and Coma sample. Plotting symbols are the same as those in Figure^ Me- 
dian error bars for each sample are shown. Dashed lines represent the com- 
pleteness limits for each sample. The solid lines are the fits to the points (see 
text). The dot-dashed line shows the relation expected for galaxies formed at 
Zf = 2 and passively-evolving to the Coma galaxies according to the models 
of Maraston (2005). 



The eastern galaxies have M/Ly consistent with passively- 
evolving stellar populations formed at Z[ orm ~ 2, and their in- 
ternal scatter is consistent with the Coma FP. They also have 
a lower fraction of SOs, 1 of 6 rather than 1 1 of 17. This can 
be explained if the western galaxies have undergone star for- 
mation at a more recent epoch than z ~ 2, while the eastern 
galaxies have yet to interact with the ICM, or have not under- 
gone a burst of star formation during their incorporation into 
the cluster. 

In BDJBC we found no evidence for recent star-formation 
in any cluster E/SOs and no spatial segregation of age- 
sensitive absorption-line diagnostics. Indeed, we found stel- 
lar populations with similar luminosity-weighted mean ages 
as those in the Coma cluster. This is inconsistent with ages 
derived from our measured M/Ly via the Maraston models. 
A complication is introduced by [a/Fe] which we found to 
be enhanced by 0.14 ± 0.03 in RX J0142.0+2131 galaxies 
and could have a systematic effect on M/Ly. If this were 
the case our data indicate that enhanced [a/Fe] decreases 
M/Ly. It has been suggested that an enhancement of a- 
element s will increase the blue luminosity of a s tellar pop- 
ulation ( Salasni ch et al. 20001 iThomas et al. 20031) . and pre- 
liminary models from Maraston (private comm.) predict Ly 
is up to 20% higher for young galaxies with solar metallicity 
and [a/Fe] = 0.3. Qualitatively, therefore, the models and 
data act in the same sense. We might expect to see this ef- 
fect as a correlation in residuals from the log M/ Ly — log M 
and the [a/Fe]— log M relations. No correlation is seen in 
either Coma or RX J0142. 0+2131 galaxies. However, the er- 
rors on [a/Fe] derived in BDJBC are relatively large (~ 1/5 
of the dynamic range), and the effect of age, metallicity, and 
internal scatter on log M/Ly are not clear. A more thorough 
analysis using a larger data set will be required to decouple 
the effects of these quantities on mass-to-light ratio. 

We speculate that the systematically lower M/Ly in 
RX J0142. 0+2131 is due to a cluster-cluster merger which 
must have occurred at z > 0.85. Rapid star formation 
episodes serve to increase [a/Fe] and so decrease M/Ly fur- 
ther. Individual galaxies interact with the cluster at different 
epochs which introduces a large scatter in the FP. The im- 
print of such an interaction should still be clearly visible in 
the large-scale distribution of gas or galaxies. A more de- 
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tailed picture can be painted with high-resolution X-ray imag- 
ing or multi-object spectroscopy of ~ 100s of galaxies in 
RX J0142.0+2131. 

If it is challenging to explain the history of stellar popula- 
tions in RX J0142. 0+2131, then it is near impossible to imag- 
ine an evolutionary path from their position at z = 0.28 to 
the Coma cluster. The [a/Fe] ratios can be decreased if there 
is a significant amount of merging between the bright E/SOs 
and galaxies which have formed their stars over longer peri- 
ods. Such mergers would have to occur in the 3 Gyr available 
without any star formation to avoid decreasing the Mj Ly fur- 
ther. 

7. CONCLUSIONS 

We have established the FP for E/SO galaxies in 
RX J0142. 0+2131. There is no evidence that the slope of this 
FP is different from that of the Coma cluster. On average, 
M/Ly for galaxies in RX J0142.0+2131 is 0.28 ± 0.03 dex 
lower than in the Coma cluster. 

The internal scatter of the FP in RX J0142.0+2131 is larger 
than for Coma. This can also be visualised as a larger scat- 
ter in mass-to-light ratio for galaxies in RX J0142.0+2131. 
Lenticular galaxies have lower average M/Ly with a larger 
scatter than Es. Spectroscopically-confirmed members east 
of 01h42m06s, including the second brightest cluster galaxy, 
have M/ Ly ratios consistent with passively-evolving galax- 
ies formed at Zf ~ 2 and an internal scatter consistent with 
the Coma FP. 

The large scatter in the FP can be explained by spatially- 
segregated galaxies with different average M/Ly, consis- 
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tent with galaxies in RX J0142. 0+2131 having undergone 
bursts of star formation over a range in epochs. This effect is 
not evident when looking at age-sensitive absorption line in- 
dices. However, galaxies in RX J0142.0+2131 have enhanced 
[a/Fe] which might serve to push any putative star formation 
episode to higher redshift. If this were the case our data sug- 
gest that increased [a/Fe] acts to reduce M/Ly. Improved 
modeling of the effect of age and [a/Fe] on M/Ly is re- 
quired before more accurate statements about star-formation 
epochs can be made. 

We speculate that the large scatter in galaxies' M/Ly is 
caused by a cluster-cluster merger. This scenario is consis- 
tent with the large cluster velocity dispersion and implies that 
gaseous or galactic substructure should be detectable with X- 
ray imaging or spectroscopy of hundreds of cluster galaxies. 
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